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AN  AB  INITIO  COMPUTATIONAL  STUDY  OF  METHANO 
AND  ETHANO  BRIDGED  DERIVATIVES  OF  OXYALLYL 


♦ 


by  Andrew  S.  Ichimura,  Paul  M.  Lahti,  and 
Albert  R.  Matlin 

Department  of  Chemistry,  University  of  Massachusetts,  Amherst,  MA  01003 

and 

Department  of  Chemistry,  Oberlin  College,  Oberlin,  OH  44074 


Abstract:  Ab  initio  SCF-CI  computations  on  methano-bridged  oxyallyl  (OA)  3  and 
ethano-bridged  OA  4  indicate  that  the  OA's  are  ground  state  singlet  species,  with  strong 
C=0  bonds  and  pronounced  diradicaloid  character  in  the  Cl  wavefunctions.  The 
computational  excited  B?  pi-pi*  states  of  3-4  are  ~4-7  kcal/mol  higher  in  energy  than 
the  ^A-j  ground  state  diradicaloids.  The  bicyclobutanone  analog  of  3S,  5,  is  predicted  to 
be  substantially  lower  in  energy  than  any  of  the  diradicaloid  forms,  but  the 
bicyclopentanone  analog  of  4S,  6,  is  somewhat  higher  in  energy  than  4S.  Experimental 
observation  of  the  triplet  manifold  of  these  OA  derivatives  is  likely  to  be  rendered 
quite  difficult  by  their  small  singlet-triplet  gaps,  and  the  thermodynamic  accessibility  of 
the  bicyclic  analogs.  Qualitative  molecular  orbital  considerations  suggest  that  electron- 
rich  olefins  might  be  trapped  by  a  nonconcerted  mechanism,  in  accord  with  recent 
experimental  results. 
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AN  AB  INITIO  COMPUTATIONAL  STUDY  OF  METHANO 
AND  ETHANO  BRIDGED  DERIVATIVES  OF  OXYALLYL 


by  Andrew  S.  Ichimura,!  Paul  M.  Lahti  }■  and 
Albert  R.  Matlin^ 

Department  of  Chemistry,  University  of  Massachusetts,  Amherst,  MA  01003 

and 

Department  of  Chemistry,  Oberlin  College,  Oberlin,  OH  44074 

Introduction:  Oxyallyl  (OA,  1)  and  related  molecules  have  been  discussed  as 
reactive  intermediates  or  transition  states  in  a  number  of  reactions,  despite  an 
apparent  dearth  of  clear  experimental  evidence  for  their  direct  observation.  In  such  a 
situation  where  experiment  has  been  less  than  clearcut,  computational  chemical 
methods  can  be  useful  to  clarify  the  problem.  Semiempirical  and  ab  initio 
computational  studies  have  been  carried  out  on  OA’s,  with  the  most  sophisticated 

published  study  to  date  predicting  a  modest  preference  for  a  ground  state  (GS)  pi,pi^ 

3  8  1 

B?  state  in  the  parent  system  1,  despite  the  possible  stabilizing  influence  on  the 

state  of  closed-shell  zwitterionic  resonance  structure  2.  The  computed  triplet-singlet 

(T-S)  gap  was  small  enough  to  leave  to  question  as  to  whether  higher  levels  of  theory 

might  reverse  the  favoring  of  a  high-spin  GS  for  OA  Tt>  addition,  consideration  of 

various  OA's  leads  one  to  query  the  effects  of  substitution  on  the  OA  state  energy 

ordering,  lest  by  apparently  innocuous  changes  in  the  OA  system  one  might 

substantially  alter  its  nature. 

[INSERT  SCHEME  1] 

Two  such  likely  OA  derivatives  of  interest  for  experimental  study  are  methano- 

bridged  OA  3  and  ethano-bridged  OA  4  and  their  ring-closed  bicyclic  forms  5  and  6. 

9-10 

Similar  derivatives  have  been  employed  in  the  study  of  trimethylenemethane 
(TMM,  7),  where  the  general  system  8  has  been  much  studied.  Since  bridging  of  7  to  8  is 
accompanied  by  strain  incorporation,  cyclization  is  less  favorable  than  for  the 


unbridged  system,  making  it  is  easier  to  study  the  diradical.  Use  of  methano  and 
ethano  bridging  in  the  OA  system  could  in  principle  lead  to  a  similar  favoring  of  the 
ring-opened  forms  3  and  4.  The  following  computational  study  is  aimed  at  elucidating 
the  GS  multiplicity  of  3-4,  the  approximate  stability  of  the  ring-open  vs.  the  bicyclic 
forms,  and  the  general  electronic  nature  of  the  bridged  OA's.  Given  recent 
theoretical^  and  experimental^  ^  interest  in  unusual  bond-deficient  species  that 
may  by  various  points  of  view  be  described  as  either  non-Kekule  or  nonclassical 
Kekule  molecules,  a  question  of  particular  interest  (also  asked  by  the  study  of  parent 
OA  1)  is  whether  these  OA  derivatives  are  well  described  as  non-Kekule  diradicals,  or 
as  zwitterions  (e.g.,  2). 

Computational  Methods:  All  computations  in  this  work  were  done  on  the 

University  of  Massachusetts  Department  of  Chemistry  Celerity  1260-D  UNIX 

computer.  Geometry  optimization  for  all  molecules  was  carried  out  using  the 

program  GAMESS.  For  the  diradicals,  only  pi-pi  occupancies  were  studied  for  open- 

* 

shell  cases,  since  n-pi  states  have  typically  been  found  to  be  substantially  higher  in 

Q 

energy  for  cases  such  as  parent  OA  1.  Unrestricted  Hartree-Fock  (UHF)  wavefunctions 
were  used  for  triplet  (T)  state  optimizations,  since  it  has  been  noted  that  they  generally 
give  a  better  account  of  correlation  effects  in  geometry  predictions  than  do  restricted 
open-shell  (ROHF)  wavefunctions.^  For  the  singlet  (S)  bicyclic  molecules  5-6, 
restricted  (RHF)  wavefunctions  were  used  for  optimization,  while  for  singlet  diradicals 
3S  -  4S  two-configuration  SCF  wavefunctions  were  used.  Computer  space 
considerations  prevented  our  use  of  multiconfiguration  SCF  optimization  with  a 
reasonable  basis  set.  For  the  diradicals,  C-2v  symmetry  was  imposed  on  the 
optimizations,  with  Cg  symmetry  used  for  the  bicyclic  molecules.  Standard  gaussian 
split-valence  basis  sets  were  used,  with  modified  polarization  functions  where  used. 
Instead  of  the  usual  polarization  d-orbital  exponent  a  =  0.80  for  all  atoms  in  the  first 


row  of  the  periodic  table,  for  optimizations  only  we  used  acarbon  =  ^  ^  anc*  aOxytTen 

1  t 

=  0.80.  In  all  cases  the  rms  gradients  were  reduced  to  one  millihartree/bohr  or  less. 

Once  appropriate  geometries  were  found  at  the  various  levels  of  theory,  single 

point  configuration  interaction  (Cl)  computations  were  carried  out  using  the  MELD 

series  of  programs,  with  standard  gaussian  basis  sets.  Excitations  up  to  quadruples 

(SDTQ-CI)  were  generated  from  the  ROHF  B2  reference  for  triplets  and  from  a  closed- 

shell  ^Aj  reference  state  that  used  the  TCSCF  MO's  for  singlets.  Instead  of  the 

canonical  virtual  orbitals,  appropriately  transformed  K-orbitals  were  used  in  the  Cl 

to  recover  correlation  energy  more  efficiently.  The  active  orbital  space  used  for 

configuration  generation  for  the  diradicals  in  most  cases  consisted  only  of  pi  symmetry 

MO's.  As  we  shall  see,  his  will  lead  to  an  important  quanititative  limitation  to  our 

results,  although  the  qualitative  findings  should  not  change.  For  the  bicyclic  cases,  a 

variety  of  mixtures  of  A’  and  A"  orbitals  was  used  for  SDTQ-CI.  The  orbital  space  used 

for  Cl  computations  is  described  in  the  tables  of  results  where  appropriate. 

3  1 

Results:  Figure  1  shows  geometric  parameters  for  the  and  A^  stated  of  3-4, 

* 

as  well  as  ring-closed  bicyclic  forms  5-6,  where  a  modified  6-31 G  basis  set  was  used  for 
optimization  as  described  earlier.  Complete  cartesian  coordinate  descriptions  of  these 
geometries  are  available  in  the  supplementary  material.  Table  1  shows  a  comparison 
of  a  selection  of  geometric  parameters  for  these  species  as  a  function  of  basis  set,  to 
allow  some  evaluation  of  the  differences  caused  by  use  of  different  basis  sets  in 
optimization. 

[INSERT  FIGURE  1,  TABLES  1  -  2  ABOUT  HERE] 

*  .  . 

Table  2  shows  standard  basis  set  6-31G  SDTQ-CI  single  point  energies  in 

hartrees  for  the  OA  diradicals.  For  B2  diradicals  3T  and  4T,  excitations  were  generated 

2  2  11 

from  the  reference  configuration  (1  b^)  {2b^)  (3b  (l^  to  ^rst  n*ne  virtuals 

and  first  six  virtuals,  making  a  total  of  10,830  spin-adapted  configurations.  For  the 

1  2  2  2 
A-j  diradicals  3S  and  4S,  the  reference  configuration  was  (1  b7)  (2b  (la^)  >  with 


excitations  to  the  first  nine  virtuals  and  the  first  six  a0  virtuals  for  a  total  of  5896 

spin-adapted  configurations.  Other  SDTQ-CI  active  spaces  were  attempted,  but  this  was 

the  largest  permitted  by  the  disk  space  resources  available  to  us.  The  potential  effects  of 

★ 

neglecting  sigma-sigma  excitations  will  be  discussed  later  in  this  paper. 

Table  3  shows  the  SDTQ-CI  energies  in  hartrees  of  the  bicyclic  ketones  5  and  6, 
using  different  Cl  active  spaces  in  each  case.  The  best  choice  of  orbital  subspace  allowed 
by  our  space  constraints  employed  the  top  three  a'  MO’s,  plus  either  one  or  two  a" 
MO's  with  a  selection  of  virtual  orbitals  (see  table).  We  feel  that  the  choices  of  Table  3 
that  use  21,235  or  40,405  spin-adapted  configurations  are  reasonable  estimates  of  the 
energies  of  the  bicyclic  species  relative  to  their  diradical  counterparts,  since  these  two 
computations  are  energetically  comparable  for  both  5  and  6.  Given  the  number  of 
generated  configurations  in  these  cases,  it  is  unlikely  that  large  energy  decreases  would 
be  realized  for  the  ring-closed  forms  by  increases  in  the  Cl  active  space. 

[INSERT  TABLE  3] 

Discussion:  The  need  for  polarization  functions  in  describing  OA's  is  evident  by 
the  difference  in  geometric  parameters  between  optimizations  with  and  without 
polarization.  Much  of  the  stabilization  of  the  ^A^  state  of  the  OA's  comes  from  a 
strong  C-O  bond,  which  requires  polarization  to  be  well  described.  An  earlier  study  of 

Q 

parent  OA  by  Borden  and  Morokuma  clearly  showed  selective  stabilization  by  use  of 

polarization  functions  of  the  ^Aj  state  with  its  short  C-O  bond.  In  our  ensuing 

discussions  we  shall  only  discuss  results  where  optimization  and  energy  SCF-CI 

* 

computations  were  both  carried  out  with  the  6-31G  split-valence  polarized  basis  set. 

*  3  1 

Comparison  of  the  6-31G  geometries  for  the  B2  and  A^  states  of  both  ethano 

3 

and  methano  bridged  OA  shows  the  greater  delocalization  of  the  B2  state,  as  expected 

8  21 

bv  analogy  to  3-21G  and  more  recent  double-zeta  with  polarization  (DZP)  MCSCF 

3 

computations  by  Borden  and  Morokuma.  The  B2  C-O  bond  lengths  for  3T  and  4T 
(123.8  and  126.2  pm  respectively)  are  quite  comparable  to  a  recently-computed  length  of 


124.3  pm  found  by  Borden  and  Morokuma  for  ^2  OA  1  at  a  DZP-CI  level  of  theory."^ 

3 

The  C-C  pi-bonds  for  the  B2  state  are  in  relatively  good  correspondence  with  the  145.3 

21 

pm  found  by  the  latter  workers. 

The  diradicals  have  very  short  C-0  bond  lengths,  118.7  and  120.0  pm  for  3S 
and  4S  respectively.  These  compare  fairly  well  to  the  122.7  pm  found  by  Borden  and 
Morokuma  for  1  at  the  DZP-CI  level.-  Our  C-C  bond  lengths  also  compare  well  to  the 

146.3  pm  found  for  OA  by  these  workers.  The  short  C-O  and  reasonably  long  C-C  pi- 
bonds  found  by  ourselves  and  the  other  workers  seem  to  argue  against  representing 
OA's  with  a  major  component  of  the  zwitterionic  resonance  structure  2,  which  one 
would  expect  to  have  shortened  C-C  bonds  and  a  fairly  long  C-O  bond.  The  presence  of 
the  short  C-O  bond  implies  localization  of  unpaired  electrons  on  the  peripheral 
carbons  of  the  OA  pi-system,  and  in  the  case  of  1  this  has  been  interpreted  as  evidence 

Q 

for  the  diradical  (rather  than  zwitterionic)  nature  of  OA. 

Comparison  of  the  four  and  five-ring  cases  shows  an  appreciably  shorter  C-O 

3  1 

bond  length  for  the  smaller  ring  in  both  B2  and  A^  states.  If  one  accepts  that  the 
geometric  results  of  Borden  and  Morokuma  may  be  directly  compared  to  ours,  in  the 
*A^  manifold  particularly  it  seems  that  r q  q  decreases  as  the  allyl  C-C-C  angle 
decreases.  This  is  logical  in  terms  of  hybridization  arguments,  whereby  a  decrease  in 
the  C-C-C  angle  will  lead  to  increased  s-character  in  the  C-O  bond  with  commensurate 
C-O  shortening. 

Energetically,  the  most  notable  qualitative  result  of  our  results  is  that  both 
species  3-4  possess  singlet  ground  states  (GS’s)  by  modest  energetic  margins.  If  one 
were  to  accept  that  OA’s  are  diradicals,  Hund's  rule  predicts  a  triplet  ground  state. 
More  sophisticated  parity-based  qualitative  models  "  also  predict  a  triplet  GS  for 
OA.  Recent  experimental  experience  for  substituted  diradicals  seems  to  indicate  that 
the  simple  splitting  of  an  orbital  near-degeneracy  by  substitution  effects  need  not  make 
a  high-spin  GS  unachievable.  Triplet  ground  states  have  been  found  for  meta- 


9 1  24:  25 

benzoquinomethane  (9),  mcfa-naphthoquinomethane^  '  (10),  and  11,“  all  oi 

which  possess  substantial  SOMO  splitting  even  at  the  simple  Huckel  level  of  theory. 
At  times,  recent  results  might  seem  to  tempt  one  to  suspect  that  any  heteroatom-for- 
methylene  substitution  in  a  diradical  will  leave  the  GS  spin  multiplicity  qualitatively 
unchanged,  as  has  been  proposed  by  at  least  one  theoretical  model  of  polyradicals. “ 

[INSERT  SCHEME  2] 

A  weakness  of  this  simple  assumption  may  be  shown  by  the  following  analysis. 
The  larger  the  energetic  splitting  between  the  SOMO’s  of  a  diradical  system,  the  larger 
the  energy  incurred  to  keep  the  electrons  spin-parallel  in  different  orbitals,  and  the 
more  favorable  the  singlet  becomes  relative  to  the  triplet.  At  some  point,  a  critical 
value  must  in  principle  be  surpassed  where  a  large  HOMO-LUMO  gap  requires  a 
singlet  GS.  The  question  is,  what  is  a  large  HOMO-LUMO  gap  and  how  may  one  relate 
this  variable  to  experiment?  In  the  area  where  the  critical  value  is  surpassed,  a  small 
difference  in  the  HOMO-LUMO  gap  changes  a  molecule  qualitively  from  triplet  to 
singlet  GS.  Experimentally,  at  the  present  we  in  almost  all  cases  can  find  with  some 
confidence  only  the  qualitative  GS  multiplicity,  not  the  quantitative  T-S  or  S-T  energy 
gap.  A  decreased  T-S  gap  in  a  heterosubstituted  diradical  such  as  meta- 
benzoquinomethane  is  expected  by  qualitative  theory,  but  cannot  experimentally  be 
probed  at  present.  It  is  thus  not  unlikely  that  a  number  of  such  species  that  would 
qualitatively  be  expected  to  have  reduced  T-S  gaps  have  still  been  observed  to  have 
high  spin  GS's.  Sooner  or  later  heterosubstituted  species  with  low  spin  GS’s  were 
bound  to  be  elucidated,  based  on  the  argument  above.  3,4-Dimethylenefuran  12  and 
3,4-dimethylenethiophene  13  are  species  that  computationally  we  feel  are  best 
defined  as  diradicals,  which  experimentally^’^  have  been  found  to  have  singlet 
ground  states.  We  propose  that,  like  these  species,  the  bridged  OA's  3-4  also  have 
crossed  the  subtle  energetic  borderland  that  separates  a  hetero-substituted  diradical 


with  a  Hund's  rule  high  spin  GS  from  a  system  with  sufficient  HOMO-LUMO  splitting 
to  have  a  singlet  GS. 

Qualitative  MO  considerations  show  why  the  bridged  OA's  are  particularly 
suited  to  possess  HOMO-LUMO  splitting  sufficient  to  favor  a  singlet  GS.  Figure  2 
shows  the  effect  of  conceptually  activating  various  orbital  effects  in  OA  relative  to  its 
hydrocarbon  analog,  TMM.  Oxygen  for  methylene  substitution  will  qualitatively 
depress  the  energy  of  the  symmetric  orbital,  leaving  the  a ^  orbital  unaffected  to  a 
crude  approximation.  Decrease  of  the  all vl  C-C-C  angle  due  to  the  constraints  of  ring 
formation  will  reinforce  this  trend  by  through-space  interactions  that  will  raise  the 
orbital  energy  and  decrease  the  b ^  orbital  energy.  Finally,  there  will  be  a  through-bond 
effect  of  the  bridging  group  in  the  OA’s,  which  is  in  opposite  directions  for  the 
methano  vs.  ethano  bridge.  This  effect  has  been  predicted  by  Schoeller  to  be  of 
potential  importance  in  the  bridged  TMM  cases  —  the  reader  is  referred  to  the 
literature*"  for  a  more  complete  description  of  the  direction  and  nature  of  qualitative 
through-bond  interactions  in  the  TMM  system.  Overall,  the  through-bond  effect 
should  be  too  small  to  overwhelm  the  combined  effects  of  substitution  and  angle 
contraction  in  the  bridged  OA's,  and  so  will  be  of  secondary  importance. 

[INSERT  FIGURE  2] 

The  previous  analysis  may  be  confirmed  in  a  crude  qualitative  manner.  We 
carried  out  RHF  singlet  computations  upon  the  singlet  GVB  geometries  of  3S  and  4S  as 
a  measure  of  HOMO-LUMO  splitting  in  the  four  vs.  five-member  ring  systems.  We 
found  at  the  6-31G  level  that  the  pi  HOMO-LUMO  gap  for  3S  was  -  (-0.302  -  -0.015)  = 
0.287  hartrees  or  7.8  eV.  For  4S  the  analogous  gap  is  -  (-0.251  -  -0.120)  =  0.131  hartrees  or 
3.6  eV.  Since  the  actual  Cl  computations  are  based  upon  TCSCF  MO's,  it  would  be 
wrong  to  overinterpret  these  RHF  numbers,  but  in  terms  of  the  crude  HOMO-LUMO 
gap  analysis  the  trend  of  increased  gap  with  decreased  C-C-C  angle  is  obvious. 


~)  1 

Borden  and  Morokuma's  results  to  date-  indicate  that  parent  O A  1  retains  a 

3 

small  preference  (-1  kcal/mol)  for  a  B7  GS,  in  agreement  with  the  qualitative 

"7 -7.-73 

connectivity  models.-'  Given  the  influence  of  angle  contraction  in  the  bridged  OA 
systems,  our  computed  S-T  gaps  of  -4-6  kcal/mol  seem  quite  in  line  with  the 
qualitative  MO  arguments  above.  Thus,  it  seems  less  surprising  that,  despite 
substantial  efforts,  there  is  to  our  knowledge  no  direct  experimental  evidence  for 
observation  of  a  triplet  OA  species.  Based  on  our  SCF  +  pi-SDTQ-CI  results,  OA  simply 
appears  to  be  a  system  where  HOMO-LUMO  splitting  reduces  the  T-S  gap  until  small 
effects  such  as  bridging  and  C-C-C  angle  may  reverse  the  "expected"  triplet  GS 
preference.  This  effect  is  reinforced  in  the  bridged  OA’s  by  the  previously  mentioned 
(vide  supra)  hybridization  effects  which  shorten  and  strengthen  the  C-O  bond, 

favoring  the  ^  A,  state  more  than  in  unbridged  1. 

* 

Sigma-sigma  excitations  were  not  used  in  our  Cl  computations  due  to  the  space 

requirements  of  such  computations  at  the  SDTQ-CI  level.  For  these  systems  the 

cumulative  effect  of  such  excitations  would  be  large.  In  computational  studies  of 

hydrocarbon  pi-diradicals,  sigma-pi  orbital  energy  separability  could  be  assumed 

28 

without  major  errors  in  relative  state  energies.  However,  a  study  of  formyloxy  14  bv 

29 

Borden,  Davidson,  and  Feller  clearly  showed  that  sigma  excitations  can  be  important 
for  even  a  qualitative  description  of  state  ordering  in  an  open  shell  system  where  a 
heteroatom  may  bear  unpaired  spin.  Such  considerations  could  confound  our  study 
with  only  pi-excitations,  due  to  the  energetically  high  lying  lone  pair  orbital  of  the  OA 
system  (Figure  2). 

[INSERT  SCHEME  3] 

Results  for  the  oxyallyl  parent  system  suggest  that,  while  one  may  question  the 

quantitative  S-T  gaps  in  our  study,  still  the  qualitative  GS  spin  multiplicity  --  the  only 

* 

property  presently  verifiable  by  experiment  --  will  be  as  we  predict.  SDTQ-CI  6-31G 
computations  by  us  on  1  with  pi-onlv  excitations  and  sigma+pi  excitations  clearly  show 


a  selective  stabilization  of  the  state  relative  to  the  ^B 2  state  for  the  sigma+pi 

21 

schemes.  Recent  higher  level  computations  on  1  confirm  this  trend.  Since  our  pi- 
onlv  Cl  results  already  predict  a  lAl  GS  for  both  3  and  4,  we  feel  that  inclusion  of  sigma 
excitations  with  their  commensurate  demands  in  CPU  time  and  disk  space  are  not 
required  to  answer  the  experimental  prediction  of  the  GS  spin  multiplicity  of  the 
bridged  OA's. 

One  may  argue  that  the  bridged  OA's  at  least  should  not  be  classified  as  diradical, 
non-Kekule  species,  but  as  some  sort  of  zwitterionic,  Kekule  species  such  as  the 
resonance  structure  2  shown  earlier.  Aside  from  the  the  evidence  of  our  optimized 
geometries,  which  seems  not  to  support  this  argument,  we  note  that  the  Cl 
wavefunctions  describing  the  ^A^  ring-opened  structures  3S  and  4S  indicate 
intermediate  diradical  character.  Borden  has  used  the  square  of  the  ratio  of  the  two 

largest  coefficients  and  C9  in  the  Cl  wavefunction  of  a  singlet  species  as  a  relative 

1*2  2 
measure  of  diradical  nature.1  Where  C^  and  C2  are  comparable  ([C^/C7]  ~  1.00),  a 

two-configuration  wavefunction  is  required  for  a  good  electronic  description  of  the 

molecule,  which  is  a  hallmark  of  a  diradical.  For  a  Kekule  species  that  is  well  described 

bv  a  single  electronic  configuration,  C^  »  C 2,  and  [C|/C2]  »  1-00.  For  the  A^  bridged 

OA's  3S  and  4S,  [C-^/C^2  is  4.16  and  4.05  respectively.  This  is  greater  than  the  1.00-2.50 

range  that  one  might  expect  for  a  nearly  "pure"  diradical,  but  still  indicates  substantial 

diradical  character.  Borden  found  that  [Q^/C^  =  1-73  for  the  putative  diradical  12,  but 

a  value  of  5.06  for  carbonyl  ylid  15. 1  By  this  comparison  our  OA's  are  more 

diradicaloid  than  carbonyl  ylid,  but  less  so  than  3,4-dimethylenefuran  and  related 

compounds. 

Alternatively,  one  can  consider  the  occupation  numbers  of  the  Cl  natural 

orbitals  (NO's),  which  can  reflect  the  diradical  nature  of  the  two  SOMO's^  without  the 

2 

extreme  sensitivity  to  small  changes  in  C|  and  C2  that  affect  the  ratio  [C^/C?]  .  We 
find  the  ^A^  ring-open  OA's  3S  and  4S  to  have  6-31  SDTQ-CI  3&|/la?  NO  occupancies 


of  1.608/0.390  and  1.549/0.401.  These  results  agree  with  the  ratio-squared  results  in 
showing  an  intermediate  range  of  diradieal  behavior  for  the  bridged  OA's.  A  "pure" 
diradical  would  have  SOMO-^/SOMC^  NO  occupation  numbers  of  1.00/1.00,  while  a 
Kekule  molecule  would  have  a  SOMO^  occupancy  of  -2.00  and  a  SOMO?  occupancy  of 
-0.00.  We  find  that  parent  OA  at  a  similar  level  of  theory  has  NO  occupancation 
numbers  in  the  ^A^  state  of  1.37/0.64,  showing  substantially  more  diradical  nature 
than  the  bridged  species.  By  our  results,  OA's  might  best  be  described  as  diradicaloid, 
but  less  so  than  a  classical  diradical  such  as  TMM. 

In  experiments  generating  diradicals,  it  has  often  been  of  importance  to  evaluate 
the  role  of  potential  ring-closed  products  that  could  deplete  the  ring-open  form 
population.  For  instance,  the  role  of  bicylopentane  8  in  the  chemistry  of  the  ethano- 

Q 

bridged  TMM  system  7  is  experimentally  quite  important.  Our  results  suggest  that  the 
ring-closed  systems  5  and  6  are  thermodynamically  quite  competitive  with  the  ring- 
open  states.  For  the  five-member  ring  case,  the  bicyclopentanone  6  lies  between  (+3.0 
kcal/mol)  the  GS  diradicaloid  A^  and  states.  For  the  four-member  ring  case,  the 
bicyclobutanone  5  lies  a  surprising  34  kcal/mol  below  the  A^  diradicaloid. 

[INSERT  FIGURE  3] 

A  study  analogous  to  ours  for  OA  was  carried  out  by  Dixon  and  coworkers  for 
32  .  33 

ethano  and  methano-bridged  TMM.  In  this  study,  the  ring-closed  8  was  found  to 

3  1 

lie  10  kcal/mol  higher  than  the  B2  diradical  TMM  7,  with  the  A^  state  19  kcal/mol 
3 

above  the  B2  (Figure  3).  By  comparison,  methylenebicyclobutane  17  was  found  to  be 
essentially  degenerate  v  ,;h  the  analogous  B2  TMM  16,  and  most  likely  lies  even 
lower  in  energy  due  to  oasis  set  level  used  in  this  computation:  the  ^A^  state  of  16 
lies  approximately  unchang  relative  to  the  B2  at  +22.0  kcal/mol. 

The  stabilization  of  the  ring-closed  form  relative  to  the  diradicals  by  methano 
vs.  ethano  bridging  for  the  TMM's  and  OA's  is  analogous,  for  reasons  that  are  probably 
based  upon  the  ring  strain  in  3  relative  to  4.  In  the  methano  bridged  system,  both  the 


OA  and  the  ring-closed  form  retain  a  substantial  amount  of  small-ring  strain  which 
will  tend  to  cancel  in  considering  the  gap  between  the  open  and  closed  forms.  In  this 
case,  the  extremely  short  C-O  bond  and  extra  C-C  bond  in  5  will  cause  it  to  be  favored 
thermodynamically  over  the  open  form  3.  In  the  ethano  bridged  system,  presence  of 
the  short  C-O  bond  and  extra  C-C  bond  in  6  will  be  offset  in  the  open  form  4  by  a  lack  of 
major  ring  strain  in  the  five-member  ring.  As  a  result,  one  might  expect  the  closed 
and  open  five-member  ring  forms  to  be  energetically  more  comparable  by  this  analysis 
than  is  the  case  for  the  four-member  ring  case,  as  is  predicted  by  our  computations. 
The  high  degree  of  strain  in  the  OA  3  and  the  relative  stability  of  the  closed  form  5 
would  seem  to  render  the  four-member  ring  OA  system  to  be  a  less  attractive 
experimental  target  than  the  five-member  ring.  In  a  similar  attempt  to  study  the 
methano-briaged  TMM  system  16,^  no  evidence  for  a  TMM  was  obtained,  but  only 
unrelated  products  apparently  derived  from  a  thermally  allowed  bicyclobutane  to 
butadiene  rearrangement  of  17. 

Zwitterionic  OA  systems  have  been  widely  invoked  as  a  key  reactive 
intermediate  in  the  photoisomerization  of  cross  conjugated  dienones  (equation  1).^ 

[INSERT  SCHEME  5] 

The  evidence  supporting  the  intermediacy  of  a  zwitterionic  species  rests  largely  on  the 
results  of  trapping  studies  with  nucleophiles  such  as  alcohols  and  halides^  The 
photoreaction  can  give  a  variety  of  deep  seated  rearrangements  and  fragmentation 
products  that  are  highly  dependent  on  both  the  structure  of  the  dienone  and  the 
reaction  medium.  An  OA  zwitterion  has  been  suggested  as  an  observable  intermediate 
that  is  stable  at  reduced  temperature  in  some  related  photochemical  solution  phase 
work.  Observation  of  the  dienone  photorearrangement  in  the  gas  phase  has  in  some 
cases  been  interpreted  as  meaning  that  the  isomerization  may  proceed  through  a 
diradical  and  not  a  zwitterionic  intermediate.^  The  zwitterionic  interpretation  of  the 
photoreaction  may  also  be  consistent  with  our  description  of  the  ground  state  of  4  as 


being  diradicaloid,  rather  than  a  true  diradical.  It  is  possible  that  a  polar  reaction 
medium  may  cause  an  altering  of  reactivity  such  that  trapping  chemistry  is  observed 
that  is  consistent  with  that  usually  interpreted  as  being  from  a  zwitterionic  molecule. 

Recently  one  of  us  has  conducted  an  olefin  trapping  study  supporting  the 

intermediacy  of  a  bridged  OA  system  (equation  2  below).  Photochemical  irradiation 

of  18  at  reduced  temperature  leads  to  a  storable  new  reactive  molecule,  believed  to  be 

the  isomer  19.  The  likely  presence  of  OA  20  is  deduced  by  the  lack  of  olefin  trapping 

products  attributable  to  trapping  by  19,  but  only  of  products  logically  derivable  from  20, 

possibly  by  a  zwitterionic  or  diradical  mechanism.  To  our  knowledge,  this  is  the  first 

example  of  the  capture  of  a  putative  OA  system  by  a  simple  mono-olefin,  although 

39 

conjugated  dienes  have  been  observed  to  give  [3+4]  adducts  with  OA  intermediates. 

[INSERT  SCHEME  6] 

The  experimental  results  of  the  olefin  and  diene  trapping  studies  is  reasonable 
in  terms  of  simple  molecular  orbital  (MO)  models  of  OA.  As  shown  earlier  in  Figure 
2,  the  OA  HOMO  will  be  of  b7  symmetry,  the  LUMO  of  symmetry.  Since  the  HOMO 
is  substantially  lowered  relative  to  the  HOMO  of  TMM,  one  might  expect  OA's  to  act 
more  as  electron  acceptors  than  TMM's,  as  observed  experimentally  (equation  2, 
above).  Equation  3  shows  the  interaction  of  the  HOMO  of  an  electron  rich  olefin  with 
the  LUMO  of  an  OA.  A  concerted  cycloaddition  cannot  occur  with  these  MO's,  hence  a 
stepwise  zwitterionic  or  diradical  path  is  likely,  as  postulated  in  experimental  studies. 
The  finding  of  fused  and  bridged  adduct- type  products  in  the  reaction  of  equation  2  is 
consistent  with  such  a  stepwise  mechanism.  As  a  result,  one  might  expect 
stereochemical  scrambling  in  trapping  studies  where  such  effects  might  be  monitored. 
Experimental  studies  are  underway  by  one  of  us  to  test  the  stereochemical  outcome  of 
olefin  trapping  by  OA's,  and  will  be  reported  in  due  course. 


[INSERT  SCHEME  7] 


Summary  and  Conclusions:  This  computational  study  of  model  bridged  oxyallyl 
derivatives  highlights  several  characteristics  of  these  species  which  are  worthy  of 
experimental  verification  as  part  of  the  present  interest  in  elucidating  factors 
determining  ground  state  multiplicity  and  electronic  nature  in  diradicals  and  related 
bond  deficient  species. 

(1)  Both  methano  and  ethano  bridged  OA's  favor  a  ground  state  with  a 
strong  C-O  bond.  Structurally,  they  do  not  appear  to  be  well-described  by  a  set  of 
zwitterionic  resonance  structures  such  as  2,  but  rather  by  diradicaloid  structure  1  with 
unpaired  electron  density  on  the  periphery  carbons  of  the  pi-system.  Efforts  to  observe 
the  low-lying  state  are  likely  to  be  impeded  by  energetically  favorable  intersystem 
crossing  to  the  singlet  manifold,  aided  by  spin-orbit  coupling  due  to  the  presence  of  the 
oxygen  atom. 

(2)  The  bridged  OA's,  although  they  possess  a  fair  degree  of  diradical  nature 
in  their  Cl  wavefunctions,  at  our  level  of  theory  seem  to  be  better  described  as 
diradicaloid  singlet  molecules,  rather  than  true  diradicals.  This  conclusion  is  based  on 
a  comparison  of  our  ab  initio  Cl  wavefunctions  to  those  computed  by  others  for 
heteroatom-containing  bond-deficient  species. 

(3)  The  bicyclic  species  5-6  are  of  potential  experimental  importance  for 
methano  and  ethano  bridged  systems.  The  stability  of  bicyclobutanone  5  is  so  great 
relative  to  the  ring-open  species  3S  and  3T  is  such  that  efforts  to  observe  the  latter  are 
likely  to  be  extremely  difficult  unless  a  high  barrier  to  ring  closure  of  3S  to  5  exists. 
Given  the  allowed  nature  of  this  cyclization,  such  a  barrier  seems  fairly  unlikely.  The 
relative  stability  of  bicyclo-pentanone  6  is  similar  to  that  of  the  ring-open  forms  of  OA 
4,  hence  the  potential  role  of  6  in  the  chemistry  of  the  five-member  ring  OA's  is  not 
clear. 
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Table  1 


Comparison  of  Selected  Geometric 

Parameters  as 

a  Function  of 

Basis  Set 

for  Bridged 

Oxyallyl  Systems 

Molecule 

3-4. 

Parameter 

3-21G3 

rx  K 

3-21  GO  ’ 

6-31 G3 

6-31G  3 

3  UHF 

C.-O 

128.2 

123.4 

128.3 

129.8 

CVC2 

cl-ci 

144.7 

144.8 

143.9 

144.0 

154.7 

152.6 

153.3 

151 .9 

C1-C2'C3 

cl-ci-ci 

92.8 

92.0 

93.0 

92.2 

85.3 

86.0 

85.8 

86.1 

°4°1  °2 

90.9 

91.0 

90.6 

90.9 

3  TCSCF 

CO 

Cf-c 

c;< 

120.6 

119.5 

121 .8 

1 18.7 

148.6 

143.7 

146.6 

146.0 

154.6 

151 .4 

153.4 

152.2 

C1C2'C3 

C2-C3-C4 

C4'C?-C2 

87.9 

85.5 

88.4 

87.8 

83.7 

80.2 

83.5 

88.4 

94.2 

97.2 

94.0 

94.4 

4  UHF 

ClO 

c‘-c 

C2‘C3 

Cn'C4 

130.9 

126.2 

130.6 

126.5 

141.8 

142.3 

141.9 

142.1 

'  151.4 

150.1 

150.8 

150.0 

157.5 

156.0 

156.5 

155.7 

C5  Cl"C2 
C1-C2'C3 
C2'C3-C4 

110.4 

109.5 

110.5 

109.7 

110.4 

110.5 

110.0 

110.3 

104.4 

104.8 

104.7 

104.8 

4  TCSCF 

CO 

CVC2 

cl-cf 

122.1 

119.1 

123.5 

120.0 

146.1 

145.9 

144.9 

144.8 

149.7 

148.6 

149.4 

148.6 

cl-c3. 

156.1 

154.6 

155.0 

154.2 

C5C1-C2 

-CJ 

103.9 

103.5 

104.5 

103.8 

94.2 

97.2 

94.0 

94.4 

83.7 

80.2 

93.5 

88.4 

aBond  lengths  in  picometers,  bond  angles  in  degrees.  Optimization  performed  as  described  in  the 


paper.  bStandard  3-21 G  basis  set  augmented  by  polarization  functions  on  carbon  and  oxygen  as  in  Harihan, 


P.  C.;  Pople,  J.  A.  Theor.  Chim.  Acta,  1973,  38,  213.. 


Table  2 


6-31G*  SDTQ-CI  Ab  Initio  Energies  of  Oxyallyl  Species  3-4. 


Structurea 

Energy^ 

Spin  Adapted 
Configurations  (Mult)c 

S-T  Gap' 

3  UHF 

-228.58415 

10,830  (T)d 

4.0 

UHF 

-228.58421 

5,896  (S) 

GVB 

-228.58185 

10,830  (T) 

GVB 

-228.59061 

5,896  (S)d 

0.0 

4  UHF 

-267.64878 

10,830  (T) 

UHF 

-267.66415 

5,896  (S)  , 

GVB 

-267.65924 

10,830  (T)d 

6.9 

GVB 

-267.67018 

5,896  (S)d 

0.0 

aSee  Figure  1  for  structural  parameters.  ^Energy  in  hartrees  (1  hartee  =  27.21  eV 
627  kcal/mol)  from  6-31G*  SDTQ-CI  computation.  cSpin  adapted  configurations 
generated  by  MELD  and  used  in  SDTQ-CI  computation  for  the  indicated 
multiplicity  at  the  indicated  fixed  geometry.  ^Lowest  energy  computed  for  the 
indicated  multiplicity.  eEnergy  in  kcal/mol  relative  to  lowest  energy  of  lowest 
state. 


Table  3 


* 


Selected  6-31G* 

SDTQ-CI  Energies  for 

ricyclic  Compounds  5  and  6. 

Spin  Adapted 
Configurations  (MO's)c 

Structure3 

Energy^ 

5 

-228.63359 

28,804  ([11a’  -  13a’]  =>  7xa') 

(  5a"  =>  4xa") 

-228.63859 

40,405  ([11a’  -  13a']  =>  7xa’) 

(  5a”  =>  5 xa") 

-228.64546 

21,235  ([11a’  -  13a'l  =>  5xa’) 

([  4a"  -  5a"]  =>  3xa") 

6 

-267.66470 

21,235  ([12a’  -  14a’]  =>  5*a') 

(  7a"  -  8a"  =>  3xa") 

-267.66537 

40,405  ([12a’  -  14a’]  =>  7xa’) 

([  8a"]  =>  5xa") 

aSee  Figure  1  for  geometries.  ^Energy  in  hartrees  (1  hartee  =  27.21  eV  =  627 
kcal/mol).  cNumber  of  spin-adapted  configurations  generated  by  MELD,  using  the 
Cl  active  space  noted.  [X  a"  -  Y  a"]  =>  Z xa"  means  the  electrons  in  orbitals  of  a” 
symmetry  numbered  X  -  Y  are  excited  up  the  the  quadruples  level  into  the  lowest 
lying  Z  a"  virtual  orbitals  (as  well  as  other  virtuals  in  the  overall  active  space). 
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CAPTIONS  AND  LEGENDS  FOR  FIGURES 


Figure  1. 

Caption:  6-31 G*  Ab  Initio  Predicted  Geometries  for  Bridged  Oxyallyls 

and  Kekule  Isomers. 

Legend:  Bond  Lengths  are  in  picometers,  bond  angles  in  degrees 

(underlined).  B_  geometries  obtained  using  UHF 
wavefunctions,  diradical  geometries  using  TCSCF 
wavefunctions,  bicyclic  geometries  using  RHF  wavefunctions. 
Full  geometric  details  are  listed  in  supplementary  material. 


Figure  2. 

Caption:  Perturbation  Effects  Upon  Frontier  Molecular  Orbitals  in 

Bridged  Oxyallyls  Relative  to  Trimethylenemethane. 

Legend:  All  effects  shown  are  qualitative  and  not  to  scale.  The 

smaller  effects  of  through  bond  coupling  by  ethano  vs. 
methano  bridging  of  oxyallyl  are  not  included.  Cf.  ref. 

27  for  a  discussion  of  this  effect. 

Figure  3. 

Caption:  Comparison  of  Relative  Energies  for  Methano  vs.  Ethano 

Bridged  Trimethylenemethane,  Oxyallyl,  and  Kekule  Isomers. 

Legend:  The  reference  zero  of  energy  in  each  comparative  scheme  is 

chosen  arbitrarily  to  allow  relative  comparisons  between 
different  species  --  no  comparison  of  absolute  energies 
between  different  species  is  implied  in  the  scheme. 

Energies  are  relative  to  the  zero  energy  in  each  species, 
kcal/mol.  bond  =  Kekule  isomer.  Energies  for  the 
trimethylenemethanes  are  from  refs  32-33. 
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